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Abstract —This work develops power control algorithms for 
energy efficiency (EE) maximization (measured in bit/joule) in 
wireless networks. Unlike previous related works, minimum-rate 
constraints are imposed and the signal-to-interference-plus-noise 
ratio takes a more general expression, which allows one to en¬ 
compass some of the most promising 5G candidate technologies. 
Both network-centric and user-centric EE maximizations are 
considered. In the network-centric scenario, the maximization 
of the global EE and the minimum EE of the network are 
performed. Unlike previous contributions, we develop centralized 
algorithms that are guaranteed to converge, with affordable 
computational complexity, to a Karush-Kuhn-Tucker point of 
the considered non-convex optimization problems. Moreover, 
closed-form feasibility conditions are derived. In the user-centric 
scenario, game theory is used to study the equilibria of the 
network and to derive convergent power control algorithms, 
which can be implemented in a fully decentralized fashion. Both 
scenarios above are studied under the assumption that single 
or multiple resource blocks are employed for data transmission. 
Numerical results assess the performance of the proposed so¬ 
lutions, analyzing the impact of minimum-rate constraints, and 
comparing the network-centric and user-centric approaches. 

Index Terms —Energy efficiency, resource allocation, power 
control, user-centric, network-centric, Nash equilibrium, rate 
constraints, QoS constraints, massive MIMO, relay networks, 
hardware impairments, 5G technologies. 


I. Introduction 

Currently, the percentage of the global world CO 2 emissions 
due to the information and communications technology (ICT) 
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is estimated to be 5% [1], While this may seem a small 
percentage, it is rapidly increasing, and the situation will 
escalate in the near future with the advent of 5G networks. It 
is anticipated that the number of connected devices will reach 
50 billions by 2020 [2], and that a lOOOx data rate increase is 
required to serve so many connected devices [3], However, it is 
also clear that obtaining the required lOOOx by simply scaling 
up the transmit powers is not possible, as it would result in an 
unmanageable energy demand, and in greenhouse gas emis¬ 
sions and electromagnetic pollution above safety thresholds. 
Instead, the data rate must be increased by a factor 1000, 
at a similar power consumption as in present networks. This 
requires a lOOOx increase of the energy efficiency (EE), i.e., 
the efficiency with which ICT systems use energy to transmit 
data [4], This is of paramount importance for operators (e.g., 
to save on electricity bills) and end-users (e.g., to prolong the 
lifetime of batteries) and thus has motivated a great interest 
in studying and designing power control strategies taking into 
account the cost of energy. 

Power control for energy efficiency can be performed in a 
centralized or decentralized manner. Both approaches are of 
interest in the context of 5G networks, for which network¬ 
centric techniques like Cloud-RAN and cooperative multi¬ 
point (CoMP) [5], as well as user-centric techniques like 
device-to-device (D2D) [6] and the use of femto-cells [7], 
have been proposed. In a network-centric approach, resource 
allocation is performed centrally and all network-nodes co¬ 
operate to maximize a common system-wide performance 
function. Instead, a user-centric approach is implemented in 
a distributed way, with the different nodes behaving in a self¬ 
organizing, and often competitive, fashion. Network-centric 
approaches typically grant better performance, with respect 
to self-organizing algorithms, but on the other hand they are 
more complex and require a larger feedback overhead to be 
implemented. An overview of recent results in these two 
directions is provided below. 

A. State-of-the-art 

A non-exhaustive list of recent works dealing with EE 
optimization in centralized networks includes [8]—[15] and 
references therein. Given the non-convex, fractional nature of 
EE, the main mathematical tool for centralized optimization 
of EE-related metrics is fractional programming - a branch of 
optimization theory that provides algorithms with polynomial 
complexity to globally maximize fractional functions with a 
concave numerator and a convex denominator [16], However, 
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even this powerful tool fails when interference-limited net¬ 
works must be optimized. This is due to the fact that the 
presence of multi-user interference makes the numerator of 
EE non-concave. A common way out to this problem is to 
rely on orthogonal or semi-orthogonal transmission schemes 
as well as on interference cancellation techniques (to fall back 
to the noise-limited case). Contributions in this direction are 
given in [9], [14], [15]. In [9], [15] multi-carrier networks 
are considered, and the global energy efficiency (GEE) of the 
system (defined as the ratio between the sum achievable rate 
and the total consumed power) is optimized using orthogonal 
or semi-orthogonal subcarrier allocation schemes. In [14], 
the authors consider a multiple-antenna system and aim to 
maximize the GEE when non-linear interference cancellation 
techniques are used. However, orthogonal interference sup¬ 
pression schemes inevitably result in a poor resource reuse 
factor and are not practical in large networks. An alternative 
approach consists in handling the interference by means of 
heuristic solutions, typically based on the use of alternating 
optimization techniques. Examples in this context are given 
by [10] in which the minimum of the individual EEs is 
maximized and also by [11], [12] where both the maximization 
of GEE and of the sum of the individual EEs are considered. 
While these approaches can operate in interference-limited 
networks, they do not guarantee convergence and/or are not 
supported by strong optimality claims. Moreover, they are 
typically tailored to the maximization of specific EE metrics. 
A first attempt to provide a unified framework to tackle EE 
optimization problems in a centralized way is given in [13], 
where polynomial-time algorithms to optimize the GEE as 
well as the sum and product of the individual EEs are provided. 

As for EE maximization through decentralized solutions, 
in [17] the authors study the Nash equilibrium problem for a 
group of players aiming at maximizing their own EE while sat¬ 
isfying power constraints in single and multi-carrier systems, 
similarly to what was done in [18] for rate maximization. A 
quasi-variational inequality approach is taken in [19], where 
power control algorithms for networks with heterogenous users 
are developed. In [20], [21] a similar problem is considered, 
with regard to relay-assisted systems, whereas single-user 
multiple-input multiple-output (MIMO) systems are consid¬ 
ered in [22]. However, all of these previous works do not 
account for rate requirements, and so the resulting users’ rates 
at the equilibrium could be fairly low. Incorporating target 
rates changes the setting drastically since any user’s admissible 
power allocation policy depends crucially on the policies of 
all other users. First results in this context are provided in [23] 
wherein Nash equilibria are found to be the fixed points of a 
water-filling best-response operator whose water level depends 
on the rate constraint and circuit power. 

B. Motivation 

The aim of this paper is to develop a unified framework for 
the analysis and design of both centralized (network-centric) 
and decentralized (user-centric) EE power allocation policies 
in a wireless network in which K transmitters (possibly) share 
N mutually orthogonal resource blocks for data transmission. 


Unlike most previous related works, we aim at maximizing 
different EE metrics while satisfying minimum rate constraints 
(or quality-of-service (QoS) requirements). Moreover, we as¬ 
sume that the signal-to-interference-plus-noise-ratio (SINR) 
experienced by transmitter k at its intended receiver on re¬ 
source block n takes the following general form: 

Ctk,nPk,n ... 

^Yk,n 2 . / , v—' (1) 

^k,n ' tyk,nPk,n A 2-^j^k ^kj,nPj,n 

where „ is the /c-th user’s transmit power over resource 
block n, whereas ak, n , are positive quantities 

that do not depend on the users’ transmit powers, but only 
on system parameters and propagation channels. In particular, 
ak, n an d (f>k,n are assumed to depend only on user k's 
channels on resource block n, while the coefficients u>kj, n 
depend on the other users’ channels on resource block n. The 
main motivation behind the adoption of (1) is that there exist 
several communication systems and technologies in which 
the SINR takes this form. 1 Interestingly, this is the case of 
some candidate technologies for 5G networks, e.g., practical 
massive MIMO networks in which the massive amount of 
deployed circuitry prevents the use of high-quality hardware 
and thus gives rise to hardware impairments [24]. The form 
in (1) arises also when imperfect channel state information 
(CSI) is available due to channel estimation errors. This is 
again a typical situation in practical massive MIMO sys¬ 
tems [25], Other relevant examples are heterogeneous, relay- 
assisted interference networks (e.g. multi-cell and/or small-cell 
relay-assisted MIMO networks [26], multi-cell and/or small¬ 
cell orthogonal frequency division multiple access (OFDMA) 
networks [20]), and relay-assisted D2D networks [27], Fi¬ 
nally, other well-established communication technologies are 
also included, such as ultra wide-band systems [28], or, 
generally, transmissions affected by inter-symbol interference 
and frequency-selective fading [29], In Section II-C we will 
describe in detail some examples of communication systems 
using candidate 5G wireless technologies, in which the SINRs 
take the form in (1). 

C. Contributions and paper outline 

The major contributions of this work are as follows: 

• A unified framework for EE optimization is developed 
for both centralized and decentralized networks with rate 
and power constraints in which the users’ SINRs take the 
more general expression in (1). This allows encompassing 
some of the emerging technologies for 5G. 

• The maximization of the GEE as well as of the mini¬ 
mum EE is considered in the network-centric case. Both 
problems are non-convex and thus hard to solve. We 
first derive closed-form feasibility conditions, and then 
exploit the tools of fractional programming and sequential 
convex optimization to develop centralized power control 
algorithms that are guaranteed to converge to a Karush- 
Kuhn-Tucker (KKT) point of the non-convex problems 
with affordable computational complexity. 

'Observe that (1) includes as a special case the SINR expression typically 
encountered in wireless communication systems, which can be obtained by 
simply letting <f>k n = 0 for all k, n. 
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• In the decentralized setting, the users in the network are 
modeled as rational, self-organizing agents that engage 
in a non-cooperative game wherein each one aims at 
maximizing its individual EE while targeting its own 
power and rate constraint. The existence and uniqueness 
of Nash equilibrium points are studied and a fully dis¬ 
tributed algorithm based on best-response dynamics is 
proposed to reach equilibrium. 

• The above scenarios are studied under the assumption 
that one or more resource blocks are employed for data 
transmission. 

The remainder of this paper is organized as follows. Sec¬ 
tion II introduces the signal model and formulates the EE 
maximization problems from both network- and user-centric 
perspectives. Moreover, Section II-C provides a detailed de¬ 
scription of some 5G candidate technologies for which the 
SINRs are expressed as in fl). The centralized and decentral¬ 
ized approaches for EE maximization in single-resource block 
transmissions are analyzed in Sections III and IV, respectively. 
The counterpart cases of multiple-resource blocks are studied 
in Sections V and VI. In Section VII the performance of the 
proposed algorithms is numerically analyzed with reference to 
case-studies inspired to 5G technologies. Concluding remarks 
are given in Section VIII. 


A. Network-centric formulation 

Based on the user-centric EE metric (2), two relevant 
network-centric performance metrics are investigated in this 
work. The network GEE ip is given by the system achievable 
sum-rate over the total power consumed in the system: 

, A EfcLl T,n=lB\og 2 {l+lk,n) 

Pc + Etii T Pfe 

with p c = ^2k- 1 Pc,k being the total circuit power dissipated 
in the network. 2 Another important energy-efficient metric is 
the minimum of the weighted EEs, defined as 

rj = min ^ w k p k . (7) 


Within the above setting, the GEE maximization problem can 
be mathematically formulated as: 


r 


max ip = max 

p£V p£V 


EfcLl En=l B lp g 2 (1 + Tfc.n) 

Pc + EfcLi lT Pfc 

( 8 ) 


whereas the weighted minimum-EE maximization problem 
can be written as: 


if = max p = max min w k r]k , (9) 

— pG*P — pG'P k=l,...,K 


II. Energy-Efficiency Problem Formulations 

Consider a wireless interference network with K transmit¬ 
ters, S receivers, and N available resource blocks (that might 
represent time or frequency bins) of bandwidth B in which 
the SINR of user k at its intended receiver takes the general 
form in (1). The EE rj k (measured in bit/Joule) of user k is 
defined as the ratio of the achievable rate over the N resource 
blocks and the total consumed power [16] 


„ A EL=i- B1 °g2( 1 + 7fe > n) 

V k — 


Pc,k + i^pic 


( 2 ) 


with p ck being the circuit power dissipated to operate 
the fc-th transmitter and its intended receiver, and p k = 
[Pfc,i>Pfc, 2 > • • • ,Pk,N] T € being the power allocation 

vector of user k over the N resource blocks. We assume that 
p k must satisfy the following (local) power constraint: 


1 Pfc -Pfc <0 


(3) 


where p k denotes user k’s maximum power. Unlike most 
previous related works, we assume that minimum achievable 
rates need to be satisfied. This amounts to setting: 


N 


5L 1 °g 2 (l+3' k,n)-0 k >0 


(4) 


n—1 


where 9 k is the target rate of user k (in bit/s/Hz/user). The 
feasible set of p/,. is thus given by: 


N 


V k =\ Pfc e R+ : l 1 p k < Pk, lo g 2 ( 1 + 7fc,n) >h\- (5) 


Accordingly, we call V = IIaLi the feasible set of p = 
[pf,p|\... ,Pa-] t e R% N . 


with {u>fc}fc non-negative weights. As discussed later, both 
(8) and (9) are non-convex fractional problems, which will 
be tackled by means of fractional programming theory 3 and 
sequential convex optimization. 

Remark 1. Observe that ip and // represent the two extreme 
points in the trade-off between global performance and fair¬ 
ness. In particular, ip can be seen as the benefit-cost ratio 
of the system, being defined as the ratio between the sum 
achievable rate and the total consumed power. However, it 
does not directly depend on the users’ EEs, and therefore it 
does not allow one to tune the individual EEs according to 
users ’ needs, as it might be useful in heterogeneous networks. 
On the other hand, maximizing the (weighted) minimum of the 
EEs allows us to achieve a fairer resource allocation policy. 
In particular, it is known that maximizing (7) yields a Pareto- 
efficient point where each quantity w k r) k is the same for all k. 
The whole Pareto-boundary can be simply achieved by varying 
the weights {u>fc}. However, this usually comes at the price of 
a performance loss in terms of benefit-cost ratio of the system. 


B. User-centric formulation 

A game-theoretic approach will be taken to solve the user¬ 
centric power allocation problem formulated as: 


max r) k 

Pk<ZV k 


max 

PkBVk 


ELl 51 °g2( 1 +7t,n) 

Pc,k + l T Pfc 


Vfc. (10) 


2 A slightly more general definition of the GEE considers the weighted sum- 
rate at the numerator. This might be useful to control the users’ individual 
rates. All the following results can be straightforwardly applied to this 
definition of GEE. 

3 For completeness, a brief background on fractional programming is 
provided in Appendix A. 
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Table I 

List of Functions and Variables 


Vk 

EE of user k 

T/c,n 

SINR in absence of interference and noise 


GEE 


Target rate of user k 

V 

Minimum-EE 

7 fe 

Minimum SINR requirement for user k 

'Yk,n 

SINR of user k over block n 

Pk 

Maximum power of user k 

Pk,n 

Equivalent channel gain of user k over block n 

Pc, Pc,k 

System and per-user hardware-dissipated power 


For later convenience, we denote by 


with 


A_ 

Pk,n — 


&k,n 

^kj,nPj ,n 


the equivalent channel gain when (pk /n = 0 , and call 


7fc,n 


A 


ftfe,n 

4*k,n 


(ID 


( 12 ) 


the SINR in the absence of interference and thermal noise (i.e., 
the maximum achievable SINR). Using the above definitions, 
( 1 ) can be rewritten as 4 


Ik nPk,nPk,n 

7fc,n n : 

7fc,n 4“ A 4 k,nPk,n 


or, equivalently. 


Pk,n 


'Yk,n (^ \ 

Pk,n \ 7 k,nj 


(13) 


(14) 


Note that 7 ^ n is a strictly increasing function of p kn as it 
easily follows observing that 


O 7 fe,n 7 k,nPk,n 

d p k, n (7 k,n+Pk,nPk,n) 2 


(15) 


C. Applications to 5G technologies 

As mentioned in the introduction, there exist several exam¬ 
ples of communication technologies in which the SINR may 
take the form in (1). Two case-studies are briefly detailed in 
the sequel. 

1) Massive MIMO: Consider the uplink (similar results can 
be obtained for the downlink) of a massive MIMO system 
composed of S cells wherein the base station (BS) of each 
cell uses M antennas to communicate with K single-antenna 
user equipments (UEs). Each UE is associated to a specific 
serving BS while interfering with all other UEs. As such, a 
double index notation is used to refer to each UE as e.g., “user 
k in cell j”. Under this convention, let us define h Uj G C M as 
the channel from UE j in cell l to BS i. Denoting by c,f. G C M 
the receive combining vector of UE k at its intended BS i, a 
lower bound (obtained using a standard bound based on the 
worst-case uncorrelated noise) of the uplink SINR of UE k in 
cell i takes the form [30] 


Pik 


7 k = 


E{c"h 


ik ll iik J 


PikVar{cRh iik j + U 


(16) 


4 Note that (as expected) when there is no self-interference (SI) (i.e., n = 

0 V7c, n) ^yk,n ^ Pk,nPk,n- 


life = cr 2 E{||c 4fc || 2 } pijE{\cf k h ll:j \' 2 }. (17) 

Assume that a maximum ratio combining (MRC) receiver is 
employed for data recovery. This amounts to setting c*fc = hu k 
where \\ r ik denotes the estimate of h, ( /, : given by 

h ilk = hiifc 4” h iik (18) 


with h uk being the estimation error statistically independent 
of h uk. Assume that huj ~ CAf(Q. d,; ; T:v/) where duj 
accounts for the corresponding large-scale channel fading and 
pathloss from UE j in cell l to BS i. If a minimum mean 
square error (MMSE)-based channel estimation scheme is used 
at the BS (with full pilot reuse) [30], then we have that 
Life ^ CA/*(0, PiifelA/) and Life ^ CJ\f (0, (dak Puk ) 
where 


Piik — 


d 


iik 


T + Li dilk 


(19) 


with t being a given parameter that depends on the pilot 
transmit power and the pilot sequence length. Under the above 
assumptions, we have that 


with 


__ Pik&k _ 

^T Pik4*k I L Plj^kj 


O-fe Piik ‘ tCfej — diljpiik 
(t>k = dukPiik + y Piik 


( 20 ) 


( 21 ) 

( 22 ) 


and < 7 j? = u 2 puk- Except for the cell index (which is only 
needed to ease understanding), the above SINR turns out to 
be in the same form of ( 1 ), with non-zero coefficients 7 fc. 

Similar results can be obtained when the system is affected 
by hardware impairments [24], [25]; for example, unavoidable 
clock drifts in local oscillators, finite-precision digital-to- 
analog converters, amplifier non-linearities, finite-order analog 
filters, and so forth. For the sake of simplicity, let us assume 
that the hardware impairments are only at the UEs. Following 
[25], the hardware impairments result in a reduction of the 
uplink signals by a factor 1 — e 2 with e being the error vector 
magnitude, and in an additive Gaussian distortion noise which 
carries the removed useful power. In these circumstances, a 
lower bound of the achievable SINR can be computed as 


Pik {1 - e 2 ) 


7 fc = 


E{c^h 


ik n Hk j 


Pik( 1 - e 2 )var{c^h iifc } +ft fc e 2 E{|c i ?h iife | 2 } + i ik 


( 23 ) 
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with i ik given in (17). Plugging c & = h-uk into the above 
equation and taking into account that in the presence of 
hardware impairments h a k ~ CJ\f(0, \/l — f' 2 pak I m ) and 
h uk ~ CAf(0, (duk ~ Vl - e 2 puk)lM), the SINR is found 
to be in the same form of (20). The above analysis can also 
be extended (as shown in the numerical results) to the case in 
which hardware imperfections are experienced at the BS. 

Remark 2. An alternative approach for evaluating the SINR 
in massive MIMO systems relies on exploiting the large 
dimensions of the network. Basically, results from random 
matrix theory are used to compute an asymptotic expression 
(hiown as deterministic equivalent) for (16) in the limit of 
M, K —t oo with £ (0,1). Such a deterministic equivalent 
turns out to be in the same form of ( 1 ) and close to the effective 
SINR even for a finite system [24], [30], 


2) Relay-assisted CoMP interference network: Consider 
the uplink of a two-hop multi-point to multi-point network 
with S BSs each equipped with M antennas and using 
N subcarriers to communicate with K single-antenna UEs 
exploiting a single-antenna amplify-and-forward (AF) relay. 
Denoting by the channel from user k to the relay on 
subcarrier n, the received signal at the relay can be written as 

4,1 = y/pi^h£lb k ,n + E VPi^filhn + tP (24) 

jAk 

where bk, n is the information symbol transmitted by UE k on 
subcarrier n and (E ~ CJ\f(0, off) is the relay thermal noise. 
The total power received at the relay on subcarrier n is thus 
given by 

^ r) =X>>Ei | 2 +^n- (25) 

i= i 

In order to avoid amplifier saturation, the received signal needs 
to be normalized by f } l n before being amplified by a factor 
y/prj, and forwarded to the receivers. The signal received at 
BS ik associated to transmitter k, over subcarrier n, is 


y k,n 


Pk,nPr,n 


,nh / k n bk,n H"~ 


Pi 


v ik,n 


(26) 


multi-cell and small-cell networks [26], as well as to relay- 
assisted D2D networks [27], 

Remark 3. Observe that in both scenarios described above, 
the expression of the receive filters c.jj ; and c k, n impacts the 
SINR expressions through the coefficients a k , 4> k , >kj an d 
the equivalent noise power of.. Therefore, any choice of the 
receive filters that does not depend on the transmit powers, 
results in a SINR expression as in (1). This means that the 
developed framework can be applied to both MRC and zero¬ 
forcing receivers whereas it does not readily apply to MMSE- 
based receivers as they depend on the transmit powers of UEs 
through the covariance matrix of the interference. 


III. Centralized Power Control in Networks with 
a Single Resource Block 


We start our analysis considering the case of a single 
resource block, which, for example, models single-carrier 
systems. When N = 1, deeper analytical insights can be 
gained compared to the case with N > 1. In particular, the 
single-resource block scenario is analytically more tractable 
and allows one to obtain necessary and sufficient feasibility 
conditions for the centralized energy-efficient optimization 
problems. On the other hand, only sufficient feasibility con¬ 
ditions can be obtained in the multi-resource block setting 
(Section V). This also applies to the distributed scenario in 
terms of more compact existence and uniqueness conditions 
of the equilibria (Section IV). Instead, the corresponding 
condition for the multi-resource block setting will be very 
cumbersome, and thus more difficult to handle. Moreover, 
the techniques to be presented in this section carry over to 
scenarios with N > 1, and are preparative for the more 
involved multiple resource block scenario. 

Setting N = 1 into (1) and (2) and neglecting the block 
index, ( 8 ) and (9) reduce to: 


r 


max = max 

p£V pOV 


T,k=i Bl °S,2 (1+7fc) 

Pc + Efc=iPfc 


(28) 


and 


where li, (: n £ C M is the channel vector from the relay to BS 
ik on subcarrier n and i k, n G C M is defined as 


* „ • £log 2 (l +7fc) 

77 = max 77 = max mm Wk - 

pov - P ev k—i,...,K p c ,k+Pk 


(29) 


tfc.n — 




Pj,nPr ,n 


p 


(r) 


u. h^b + 

“ifci n n j,n u j,n T 


Pr. 


P, 


(r) h ik ,nCk r) (27) 


whereas w ifci „ ~ CNf(0, of k „Im) is the thermal noise at 
receiver ik on subcarrier n. After linear reception by the filter 
c k, n and upon plugging (25) into (26) the SINR takes the form 
in ( 1 ) with 


wherein p = [pi,p 2 , ■ ■ ■ ,Pk] T G R+, 

_ _ a-kPk _ 

a k + 0 kPk + UkjPj 

and V = JlfcLi Pk with 


Vk = {Pk € K+ : Pk < P k , log 2 (l + 7k) > ^ fc }- 


(30) 


(31) 


&k,n — Pr,n\h’k,n\ I ^-'k,r^'i,k ,n I , fik,n — ll^fc.nll 

= ^Pr,ra|Cfc,nhjfc.nl + &i k || c fe,n || 

and °i,n = ^l(Pr,n\c^ n h iktn \ 2 + of k J\c k . n \\ 2 ). We should 
stress that the above results apply to any interference network 
in which the transmitters reach the receivers via an AF relay. 
For example, the above scenario applies to relay-assisted 


For later convenience, we define 

lk = 2- fe — 1 (32) 

the minimum SINR requirement for user k. Observe that 7 ^ 
must be such that 

0 < 7 fc < 7fc Vfc. 


(33) 
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The above condition follows observing that when the noise is 
negligible (i.e., —► 0 ) and only transmitter k is active then 
(30) reduces to y k = a k / 4>k — 7 k and thus the rate constraint 
log 2 (l + 7 k) > 0_ k can be met only if (33) holds true. 


A. Feasibility 

The feasibility of (28) and (29) simply amounts to verifying 
that for given values of {a k } k , {<Pk}k, and {w k ,j}k^j, the 
feasible set V is not empty. Closed-form necessary and 
sufficient conditions for ‘P to be non-empty are provided in 
the following result. 

Lemma 1. Let F £ C AxK be a matrix whose ( k,j)-th 
element is defined as 


[F] 


A_ 

k,j ~ 


0 

OCk—<t>kl_ k 


j = k 
j ± k 


(34) 


and denote by pp its spectral radius. The solutions to (28) 
and (29) exist if and only if 


Pf < 1 and (I — F) 1 s < p (35) 


where p = [p 1 ,p 2 , ■ ■ ■ ,Pk] T € M+ Xl and s £ has 
elements given by [s] fc = 0^7 k (a k - 

Proof: Following the same steps of the proof of Lemma 2 
in [31] allows us to prove the sufficiency of (35). For the 
necessity, assume that there is a vector p' satisfying the target 
rates but such that pi- > pj for some j. Since y k is a stricly 
increasing function of p k , we have that 

a kPk 

2 1 i 1 v- / + 7 k (36) 

+ fikPk + Ej^fc ukjPj 

for any p k < p k . Otherwise stated, there exists no power p k < 
p k such that 7fc = 7fc • Now, suppose that pp > 1. Since F 
is non-negative, from [32, Theorem 2.1] it follows that there 
does not exist a power vector p > 0 such that j k = fi k . This 
proves that (35) is also necessary. ■ 


B. GEE maximization 

As described in Appendix A, fractional programming pro¬ 
vides efficient tools to maximize ratios in which the numerator 
is a concave function, the denominator is a convex function, 
and the constraint set is convex, whereas no low-complexity 
optimization method is available if any of these properties is 
not met. Unfortunately, the objective function in (28) does not 
have a concave numerator, and therefore finding the global so¬ 
lution of (28) with affordable complexity appears difficult. To 
overcome this difficulty, we integrate fractional programming 
theory with the framework of sequential convex programming 
[33], This allows us to develop a computationally-efficient 
algorithm which is guaranteed to converge to a first-order 
optimal solution of (28). The general idea of sequential convex 
programming is to find local optima of a difficult problem 
with objective / to maximize, by solving a sequence of easier 
problems with objectives {/)},;. In the generic i-th step of the 
sequence, we require the following three properties: 

1 ) /;( x ) < /(x), for all x; 


2 ) fi^-V) = /(x^-D); 

3) V/iCx^- 1 )) = V/(x(‘" 1 )). 

wherein x^ _ l 1 denotes the maximizer of f, _ 1 . This approach 
has been used for resource allocation in wireless networks in 
[34], [35] for rate maximization, and, more recently, in [13], 
[26], [36] for EE maximization in different settings. 

The critical issue of this approach is to find suitable approxi¬ 
mations {fi}i which fulfill the listed requirements, while at the 
same time resulting in simpler optimization problems. As far 
as GEE maximization is concerned, this can be accomplished 
by leveraging the following lower-bound of the logarithmic 
function [37]. Specifically, V 7,7 > 0 we have that 

log 2 (l+ 7 ) > alog 2 7 + 6 (37) 


with 


a = 


7 


b = log 2 (l+ 7 )- 


7 


log 2 7- 


(38) 


1+7 1+7 

The right-hand side (RHS) and left-hand side (LHS) of (37) 
are equal at 7 = 7 , and the same holds for their derivatives 
with respect to 7 evaluated at 7 = 7 . Therefore, we may 
lower-bound + as follows: 

K K 

y, B [o fc log 2 (7 k) + b k \ y B [b k + a k log 2 ( a k pk )] 




k =1 


k =1 


Pc 



Pc 


K 


k =1 


afclog 2 ( al+fkPk + yuJkjPj 
jAk 


Pc 


K 

y?k 

&=1 

from which, letting pk = 2 qk , one gets 


K 

X> 

*;=1 


- ft (39) 


Ef=i B i b k + a k log 2 (a k ) + a k qk] 


= 


Pc + Ef=i 2« 


Ef=i B 


a k log 2 (a k +fik2 qk +E j^k 


■ (40) 

Pc + Efc =1 29fc 

Using the above results, the solution to (28) can be lower 
bounded as 


ip* >ijj* = max tp (41) 

qeQ 

with q =[(j l7 g 2 ,..., q K ] T , Q = \[ k=1 Q k and 

Qk = {qk £ R : 2 qk < P k , l°g 2 (l + 7 k) > Q_ k }. (42) 

Observe now that for any given {afc}fc and {b k }k, the nu¬ 
merator and denominator of (40) are both differentiable, and 
respectively concave 5 and convex in Finally, the set 

Q k can be shown to be convex for all k. Indeed, the fc-th rate 
constraint can be equivalently rewritten as 


2 qk {a k -f k <t>k)> 7 


—k 



( 43 ) 


5 Recall that the log-sum-exp function is convex [38]. 
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Algorithm 1 Network-centric EE maximization for N = 1 
1 : Test feasibility by Lemma 1. 

2: if Feasible then 

3: Set i = 0 and choose any p ,(J * £ “P; 

4: Set 7 ® = 7 fc°^(p < ' 0 ' ) ) and compute a®, 6 $. 0) as in (38); 

5: repeat 

6: i = * + 1; 

7: if GEE then 

8 : Solve (41) with parameters aj7 _1) and bj . 1-11 and set 

= argmax^i, p = 2 9 fc \ 

9: end if 

10: if Minimum EE then 

11 : Solve (45) with parameters a^ and b^ 1 and set 

{qi l) }k = argmax?)., p^ = 2 9 <= > ; 

12: end if 

13: Set 7 ^ = 7fc(P ( ^) and compute a^\ as in (38); 

14: until convergence 

15: end if 


Since (33) must hold true, one gets (applying the logarithm 
function to both sides) 

q k - log 2 |crfc + ^ j +log 2 j (44 ^ 

which turns out to be a convex constraint. As a consequence, 
(41) is a fractional problem which can be globally and 
efficiently solved by means of fractional programming tools 
[16], such as the Dinkelbach’s algorithm [39]. This leads to the 
general iterative procedure formulated in Algorithm 1 whose 
convergence is proved in Appendix B. 

Proposition 1. Algorithm 1 monotonically increases the GEE 
value and converges to a point fulfilling the KKT conditions 
of the original problem (28). 

C. Weighted Minimum EE Maximization 

The key difference between (28) and (29) is that the 
objective function 77 in (29) involves K fractional functions 
{rjk} rather than a single one. This makes (29) fall within 
the framework of generalized fractional programming, which 
studies the maximization of functions of multiple ratios. In 
this more general scenario, Dinkelbach’s algorithm fails, even 
assuming that each ratio {p k } has a concave numerator and a 
convex denominator. Instead, the problem can be tackled using 
an extension of Dinkelbach’s algorithm known as Generalized 
Dinkelbach’s algorithm (see Appendix A), which is guaranteed 
to converge to the global solution of a max-min fractional 
problem with limited complexity, provided each ratio has a 
concave numerator and a convex denominator [40]. Next, we 
show how the generalized Dinkelbach’s procedure together 
with sequential convex optimization can be successfully ap¬ 
plied to solve (29). 

To begin with, observe that the min(-) function is increasing 
so that the inequality in (37) can be used to lower-bound the 
solution to (29) as 

V* > 7?* = max min w k rj k (45) 

— — q GQ k=l,...,K 


where 


rj k = B 


[b k + a k log 2 (a k ) + a k q k } 


Pc,k 


- 


a k 


log 2 (c 


■ 2i k 


Pc,k 


2^ k 


(46) 


and q k is still given by q k = log 2 p k . Since each ratio in 
(46) has a concave numerator and a convex denominator, fj* 
can be computed by means of the Generalized Dinkelbach’s 
algorithm, and the maximization of 77 can be tackled as in 
Algorithm 1, whose convergence is stated in the following 
proposition and proved in Appendix C. 


Proposition 2. Algorithm 1 monotonically increases the value 
of t) and converges to a point fulfilling the KKT conditions of 
the epigraph-form representation of the original problem (29). 


Remark 4. Algorithm 1 can be straightforwardly specialized 
to maximize the system sum rate and the minimum of the users’ 
rates, since these two metrics coincide with the numerator of 
the GEE and with the minimum of the numerators of the users’ 
EEs, respectively. 


IV. Distributed Power Control in Networks with 
a Single Resource Block 

A decentralized power control algorithm looks for the 
solution of the following coupled problems [17], [20], [23]: 

arg max p k (p k ,p^ k ) Vfc (47) 

where p_ fc = [pi, • • • ,Pk-i,Pk+i, ■ ■ ■ ,Pk] T is the inter¬ 
ference vector containing all powers except user k’s, and 
V k (p- k ) is defined as in (31). This problem can be formulated 
as the non-cooperative game in normal form; 

Q = {IC,{P k (p- k )} k ,{rj k (p k ,p- k )} k } (48) 

where (in game theory parlance) K, = [1,2 ,...,AT] is the 
set of players, V k (p- k ) is player k’s strategy set, u k ( p) = 
? 7 fc(p fc ,p-fc) is player k’s utility function. The K coupled 
problems in (47) define the best-response dynamics (BRD) 
of the game, while the solution of the fc-th problem in (47) is 
the /c-th player’s best-response to the other players’ choices. 
More formally, let us define the best response B k (p-fc) of 
player k to an interference vector p_/ i; (or, equivalently, fj k as 
easily follows from (11)) as 

B k (p-fc) = arg max q k ( p k , p_ fe ) . (49) 

Pk&V k ( p^k) 

Any fixed point of the BRD is a Nash equilibrium of the game. 
In general a non-cooperative game might admit zero, one, or 
more equilibria, and even if one or more equilibria exist, the 
convergence of the BRD is not guaranteed. As a consequence, 
crucial issues in the analysis of a non-cooperative game are 
to establish the existence and uniqueness of an equilibrium, 
and whether implementing the BRD eventually yields an 
equilibrium. In our scenario, answering these questions is 
more challenging due to the fact that, unlike regular non- 
cooperative games, not only the utility functions, but also 






the players’ strategy sets are mutually coupled, depending 
on the other players’ actions p _ k . A similar non-cooperative 
game is termed a generalized non-cooperative game, and 
more restrictive conditions have to be fulfilled for a (unique) 
generalized Nash equilibrium (GNE) to exist and for the BRD 
to converge. To begin with, the following result is given: 


Lemma 2. If 


Pk>l k 


ukjPj 

a* - (t>kj_ k 


Vfc 


then B k (p~k) takes the form 

B k (p-k) = min jp fc , max j p* k ,p k j j 


wherein 


P fc (P-fc) = 


and 


A Ik 

Pk 


1 - 


Ik 

Ik 


-1 


= lk 


Oik ~ 4>klu 


Pk ~ ar § max Vk (Pk, P-fc) • 

Pf=6R+ 


(50) 


(51) 


(52) 


(53) 


Proof: The first part of the thesis easily follows from 
rewriting the rate constraints 7 k > 7 fc (using (30)) as 


Pk > l k 


+ E &k w kjPj 


Oik - fkl k 


Since p k < p k for all k £ JC, then 

+ E^fc u k jPj o\ 
-— Tk~ 


- k a k - f k 7 fc 


' Ej^fc UkjPj 


Ok - 0fc7 fc 


(54) 


(55) 


Hence, if Vfc € /C (50) holds, then there always exists a power 
p k £ [0,p fc ] such that y k > 7 fc is fulfilled. The last part of the 
proof follows by leveraging [ 20 ], where it is shown that for any 
given p /// i; is unimodal and thus admits a unique maximizer 
p k £ R+. Accounting for the power and rate constraints and 
imposing (50) eventually yields (51). ■ 


A. Analysis of the Equilibria 

The existence and uniqueness of the GNE points of Q are 
now studied under the assumption that (50) holds. 

Proposition 3. The game Q admits a nonempty set of GNE 
points. 

Proof: Observe that the existence of a GNE is guaranteed 
under the following assumptions [41]: 

1) The players’ feasible action sets Pk(P-k) are nonempty, 
closed, convex, and contained in some compact set C k 
for all p _ fc e V- k = Yl^kVe- 

2) The sets Vk{p-k) vary continuously with p_^ (in the 
sense that the graph of the set-valued correspondence 
P-fc "Pfc(p-fc) is closed). 

3) Each user’s payoff function pk{Pk, P-fc) is quasi-concave 
in p k for all p_fc £ V-k- 

In our setting, if the sufficient condition (50) is satisfied, then 
the sets Vk(p~k) are nonempty, convex , 6 closed and bounded 

6 Note that the constraint function log 2 (l + "fk) is concave in p k . 


for every p-fc. Moreover, each of them varies continuously 
with p_fc since the rate constraint log 2 (l + 7 fc) > 9 k in 
■Pfc(p-fc) is itself continuous in p-fc. Finally, following [20] 
Vk(pkiP-k) is proved to be strictly pseudo-concave since it 
is given by the ratio between a strictly concave and an affine 
function. Since any strictly pseudo-concave function is also 
quasi-concave [16], the third condition is fulfilled. ■ 

The following result shows that a unique GNE exists, and 
that the BRD always converges to such point. 

Proposition 4. The game Q admits a unique GNE point, 
which can be obtained by starting from any feasible power 
vector {pk\k=i au d iteratively updating the transmit powers 
according to (51). 

Proof: The proof builds upon the standard function frame¬ 
work [42], which states that a non-cooperative game admits 
a unique equilibrium (reachable by iteratively computing the 
players’ best-responses) provided the game admits at least 
one equilibrium and the best-response function is a standard 
function . 7 Since we have already shown that the game admits 
a GNE (see Proposition 3), we are left with proving that (51) 
is a standard function for all k. Towards this end, p k (p~k) is 
proved to be standard in [20, Appendix A], As for the function 
p k (jp_k) in (52), it is non-negative because 7 fc < j k , and it 
also fulfills the monotonicity property because it is increasing 
in all {Pj}jjtk- As for the scalability property, take any 13 > 1, 
then it holds 


P k (/3p-k) = f3f k 


+ Yljjtk^kjPj 

Ok - <t>kf k 


<Plk 


+ Ej^fc UkjPj 

o k - 4>kj k 


(56) 

%.(P-fc) • (57) 


Finally, since both p k (p_ k ) and p ( p_fc) are standard func¬ 
tions, and since p k does not depend on p we may conclude 
that (51) is also a standard function because both max(-) and 
min(-) are increasing functions. ■ 


B. Distributed implementation 

The best response of a generic player k is characterized in 
the sequel to come up with an iterative algorithm that allows 
each player to reach the GNE in a distributed manner. Toward 
this end, let us first define 


vk {x) = 7 fe 


i + 


2 Bp k 


(7fc - 9k {x)) 


(58) 


and 


9k (x) = 



(1 + 7 fc) 


(59) 


with fik and given by ( 11 ) and ( 12 ), respectively. 


7 Recall that a vector function g(p) is standard if it fulfills the properties 
of i ) non-negativity: g(p) > 0 for all p; ii) monotonicity: g(pi) > <?(P 2 ) 
for all pi y P 2 ; Hi) scalability: g(/3 p) < f3g( p), for all p and (3 > 1. 
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Algorithm 2 Iterative algorithm to solve (47). 

1 : initialize i = 0 and Vfc pk [0] € R+ in the feasible set 

2 : repeat 

3: for fc = 1 to K do 

4: receive 7 *. [i] from the serving access point 

5: compute /Zfc [*] using (64) 

6 : use fik [i] to update [t] in (63) 

7: use in (61) to run the Dinkelbach’s algorithm 

8 : set A]5[i] equal to the Dinkelbach’s output and update the 

power as: 

Pk[i + 1] = min|j5 fc ,max|7r fc j 

9: end for 

10: update i — i + 1 

11 : until convergence 


V. Centralized power control in networks with 

MULTIPLE RESOURCE BLOCKS 
In this section, we turn our attention to the case in which 
each transmitter can use multiple resource blocks, i.e., N > 1. 
Differently from the case in which N = 1, the rate constraints 
in (8) and (9) are not in a convex form. Nevertheless, the 
methodology used in Section III can be successfully extended 
to find sufficient feasibility conditions and to derive low 
complexity algorithms that converge to KKT points. 

A. GEE maximization 

Following the same steps of Section III-B, we leverage the 
lower bound in (37) to obtain 

K N 


Lemma 3. For any given p^ k (or, equivalently, p k ), the 
solution to (53) is found to be 


Pk = Ttk (Afc) = 


yk (K) 

Pk 


1 - 


Vk {X% 
7 k 


-1 


(60) 


where \* k is obtained through the Dinkelbach’s algorithm as 
the solution of the following equation: 


ip > ip 


+ & k,n lo§2 “1“ Q>k,nQk,n\ 

fc=1 n= 1 _ 

K N 

pc + E E 29k ’ n 

k—l n— 1 


K N 

( Y 

EE* 

^k,n lo§2 I ^k,n ~^~4 > k,n^ k,n “1“ ^ ^ ^kj,n^^' n 1 

k—l n—1 

V /_ 


K N 


pc + E E 


B^og 2 (1 + Vk (Afc)) — A k ( p c ,k + (A k )) — 0. (61) 

Proof: The proof is given in Appendix D. ■ 

Denote by p k [1] the transmit power of the k- th player at the 
i-th iteration step. By virtue of Proposition 4 and Lemma 3, 
it follows that an iterative algorithm operating according to 

p k [i + 1] = min|p fc ,max|7r fc (A£[*]) ,p fc [i]}} (62) 

where P k [i] is computed as (using (52)) 


k—l n— 1 

with q ktn = log 2 Pk,n- Although the numerator and denom¬ 
inator of in (65) are again jointly concave and convex in 
{qk, n }k,n (as in the case of a single resource block), the rate 
constraints in (8) are not in a convex form yet, due to the sum 
over the multiple resource blocks. To overcome this issue, we 
resort to the same ttick and lower-bound the LHS of the rate 
constraint by the concave function. 


N 

E 

n—1 


log 2 (l 


N 

' c)k,n) f. 'P \ [hk,r. 

n— 1 


p k m = 


—k 
Pk [*] 


(i-l -1 

N 

( 63 ) YE 

ak,n log 2 j 

[ <T k,n^~ < Pk,n2 qk ' n + y ]<^kj,n‘2 qj rl j 

V ik) 

n—1 


\ ) _ 


' a k,n (^§2 (t^k,n) T tlk,nqk,ri 


= Rk- 


converges to the unique GNE of Q , with p k [i\ being the 
equivalent channel gain in (11) at the i-th iteration step. The 
pseudo-code is reported in Algorithm 2. 

A close inspection of (58) - (61) and (63) reveals that the 
computation of p k [i + l] through (62) only requires knowledge 
of pk[i\- Although not available at the fc-th terminal, this 
information can be easily acquired taking into account that: 


Pk[i] 


7 k[i\ L _ 1 

Pk[i] V 7 k ) 


(64) 


where [*] denotes the SINR of Uansmitter k measured at 
its intended receiver at iteration i. Since p k [i] and f k are 
locally available at the transmitter, the computation of pk[i\ 
only requires knowledge of y k [)]. The latter can be easily 
estimated at the receiver and sent back to the corresponding 
transmitter via a return downlink channel. Therefore, besides 
being guaranteed to converge to the unique GNE, Algorithm 2 
can also be implemented in a fully decentralized fashion. 


Finally, we may write 


ip* > i/j* = max ip 

q62 


with Q — nf = i Qk an d Qk being now given by 


N 


Si- = qt e R" : ^ 2 9fc '" < p k , R k > 


Vk 


n —1 


( 66 ) 

(67) 

( 68 ) 


The solution to the above problem can be computed by means 
of Dinkelbach’s algorithm, which leads to the power allocation 
procedure illustrated in Algorithm 3. Using similar arguments 
as in Proposition 1, the following result can be proved: 8 


Proposition 5. Algorithm 3 monotonically increases the GEE 
value and converges to a point fulfilling the KKT conditions 
of the original non-convex problem ( 8 ). 


8 As observed in the proof of Proposition 2, the sequential convex optimiza¬ 
tion tool allows us to find a KKT point of the original problem also when the 
constraint functions are lower-bounded together with the objective. 
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Algorithm 3 Network-centric EE maximization for N > 1 


1: Set i = 0 and choose a starting point p (tl) ; 

2: For any k and n, set = 7fc°,l(P ( °*) and compute a} 0 } and 
b^ n as in (38); 

3: if (67) or (71) is unfeasible then 
4: break 

5: end if 
6: repeat 
7: i = i + 1; 

8: if GEE then 


Solve (67) with parameters a} 1 n 1 


anu tv- 





tit} 

= arg max p^ n = 

2 ^/n; 


10 

end if 




11 

if Minimum EE then 



12 

Solve 

(71) with parameters 

a t and 

b t and set 


tit} 

= arg max 77 ., p^ n = 

2 q t ; 


13 

end if 




14 

3 ^ 

II 

?r 

ra (p 6 )) and compute a 

t and & fe, 

\ as in (38); 


15 


B. Weighted minimum EE maximization 

A similar approach as in Section V-A can be used to solve 
(9). In particular, exploiting the fact that the min(-) function 
is increasing, using (37), and setting p kn = 2 qk - n allows us 
to lower-bound rj as 


7) > min w k f/ k =fi 
- k=l,...,K ~ 


(69) 


where 


E)„ — i B [b k n -T (X k , n log2 (&k,n) V tl k ^ n Q kjl 
IN 


Vk = 


Pc,k + E„=i 2qk ’ n 

T,n= 1 B a k,n log 2 (c^l n + 4>k,n^ k ' n + E & k W kj 


Pc,k + En=l 2<Zfc 


(70) 


Then, we have that 


rf > ft* = max min w k fi k . (71) 

— — q GQ k=l,...,K 

Observe that fj k in (70) has a concave numerator and a 
convex denominator, meaning that (71) can be globally solved 
by means of the Generalized Dinkelbach’s algorithm. The 
resulting power allocation procedure is given in Algorithm 3. 
By a similar reasoning as in Proposition 2, the following result 
can be proved. 

Proposition 6. Algorithm 3 monotonically increases the t) 
value and converges to a point fulfilling the KKT conditions 
of the epigraph-form representation of the original noil-convex 
problem (9). 

Remark 5 (Feasibility of ( 8 ) and (9)). Observe that Line 3 
in Algorithm 3 is also a sufficient feasibility test for ( 8 ) and 
(9) since both are guaranteed to be feasible provided (67) 
and (71) are feasible, which can be checked by means of a 
convex feasibility test. Moreover, this implies that (67) and 
(71) will remain feasible for all iterations of Algorithm 3. 
Obsen’e also that similar necessary and sufficient feasibility 


conditions as in Section III can in principle be derived if per- 
resource-block QoS constraints are considered. This amounts 
to enforcing log 2 (l + "/ k ,n) > 9 k ,n with 0 k , n being such that 
E^Li ®k,n — @k> f or k, with { 0 fc}fc given target rates. 


VI. Distributed Power Control in Networks with 
Multiple Resource Blocks 


As done for the single resource block case, we define 


k'k.n ( X) = 




and 


2 Pk,n + -g (7 k,n ~ 9k,n (x)) 


gk,n (x) = \/7 k, n + 4Bfik ’ n (1 + Ak,n) • ( 7 3) 


n + 


(72) 


Let us also define p ; as the power vector minimizing the 
transmit power while satisfying the rate constraints. Mathe¬ 
matically, we have that: 

p. = arg max l T p fc (74) 

p fc 6K" 

N 

subject to Y l°g 2( 1 + 7 k,n) - Ik > 0 

n—1 


from which (using the same arguments of Appendix D) one 
gets: 


Pk,n = nk ’ n Ufc) Vn 

with X k being such that: 

N 

Y log 2 (1 + v k}U (A fc )) -6 k = 0 


n—1 


and 


'Kk,n {pt) — 


Vk,n (X ) ( V k , n (x) 


l^k,r] 


1 - 


7 kf 


-1 


(75) 


(76) 


(77) 


Lemma 4. For any given p- k , the entries of 13 k (p- k ) G 
are found to be 

[B k { p_ fc )] n =min|p fcn ,max|7r fc!n (A^),p fcn || (78) 


where X k is obtained through the Dinkelbach method as the 
solution of the following equation: 


N 


N 


B\og 2 (1 + v k , n (A^))-A^ p Cjfc 


, Tl/c.n (Afc) 


n— 1 


n—1 


= 0 . 

(79) 


Proof: The proof relies on similar arguments (omitted for 
space limitations) of those in Appendix D for N = 1. ■ 

Observe that (78) can be equivalently rewritten as (since 
TTfc.n (•) in (77) is a strictly decreasing function): 

[Bfc(p -k)] n = 7 T k ,n (Afc) (80) 

with 


A^ = max {Afc, min {A^, Afc}} (81) 

and Afc such that \/n n k ^ n (\ k ) = p k n . A sufficient condition 
for (80) (or, equivalently, for (78)) to be a contraction mapping 
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Algorithm 4 Algorithm to reach the GNE of Q for TV > 1. 

1: initialize i = 0 and \/k,n p k , n [0 ] £ R+ in the feasible set 

2 : repeat 

3: for k = 1 to K do 

4: receive { 7 &,n[*]}n from the intended receiver 

5: compute {nk,n[i]}n using (84) 

6 : use {fj,k,n[i]}n to compute X k in (76) by means of the 

Dinkelbach algorithm 

7: use \ k to update i n (75) 

8 : use {fj,k,n[i]}n in (79) to run the Dinkelbach algorithm 

9: set \ k [i] equal to the Dinkelbach output and update the 

power as: 

Pk,n[i+ 1] = min max {?!>,„ (M[i\) ,p k „[*]}} 

10: end for 

11 : update i = i + 1 

12: until Vfc, n Pk,n[i] =Pk,n[i - 1] 


is provided in the following proposition (see Appendix E for 
the proof), thus ensuring the existence and uniqueness of a 
GNE for Q when TV > 1, and the convergence of the BRD. 


Proposition 7. The EE game Q admits a unique GNE and the 
BRD converges to such equilibrium whenever, for all k 


K N 

Y Y U kj,n 

7 — 1 71= 1 


sup 

I 


E 

e,n&St 


7 k,: 


9k,r. 


-£k,e 


< 1 (82) 


where El k = ( 0 ,a fc>n ] and £ k ,e = ^7 whereas 

S k ~ { n = 1) • • ■ i N : p k ,n > AJ.} and s k , n is defined as 

(Afc) — (2 + 7fe,n) 7 k,n9k,n 


A — 9k, 

s k,n — 7 fc, 


2 OLk,n (l + 7 k,n) ^' k a k,n9k,n (Afc) 


■(83) 


Similarly to the single resource block case, denote by p k , n [i] 
the transmit power of the fc-th player over block n at the t-th 
iteration step and, accordingly, compute p k , n [i\ as follows: 

-1 




Pk,n[i\ 


1 - 


lk,n[l\ 

7 fc 


(84) 


with 7 fc, n [*] being the SINR of transmitter k over block n 
measured at its intended receiver at iteration i. From the 
results of Proposition 7, it follows that the iterative procedure 
illustrated in Algorithm 4 converges to the unique GNE of Q, 
and can be implemented in a fully decentralized fashion. 



P [dBW] 


Figure 1. K = 5; M = 50; r = 10 2 . Probability of feasibility Vf versus 
P with minimum per user-rate constraints: (a) R = 15%; (b) R = 20%; 
(c) R = 25%; (d) R = 30%. 


hardware impairments are present at the UEs. However, we 
assume that hardware impairments are present at the BS, and 
we denote by cbs the resulting error magnitude. As observed 
in Section II-C1, this scenario results again in an SINR as in 
(1). Indeed, following [43] and by similar steps as in Section 
II-C1, the T'-th UE’s SINR is written as in (1) with 

a k = |hf h fc | 2 , cf) k = c 2 bs hfD fc h fc (85) 

Ukj = h"h, | 2 + e 2 BS hf Djh fc (86) 

with Dj = diag({|hj(m)| 2 }£f =1 ) and a\ = cr 2 hfh fc . 
All channels are generated according to the Rayleigh fading 
model with path-loss model as in [44]. All UEs have the 
same maximum feasible power p k = P Mk and hardware- 
dissipated power = 10 dBm Vfc. The receive noise power 
is cr 2 = FBAfo, with F = 3dB and TVo = —174 dBm/Hz. 
The minimum rate constraint 0 k is set as a percentage R k % 
of the maximum rate that user k can achieve when p k —> 00, 
while the other users’ powers are finite, namely: 

e ‘ = ® lofe(1+7 * ) = ®i log2 ( 1 + S) (87) 


VII. Numerical Results 

Numerical results are now given to assess the performance 
of the proposed solutions. To this end, two case-studies 
are considered, namely, a hardware-impaired massive MIMO 
system and a multi-carrier relay-assisted interference network. 

A. Hardware-Impaired Massive MIMO System 

Consider the uplink of a massive MIMO system as described 
in Section II-C1, with K = 5, S = 1, and M = 50. The 
communication bandwidth is B = 1 MHz, and MRC is used 
for data recovery under the assumption of perfect channel 
estimation. We also set e = 1, thereby considering that no 


For simplicity, we assume R k % = R% for all users k. 

We begin by analyzing the feasibility probability Vf of the 
EE maximization problems as a function of P for different 
values of R. The results are obtained by averaging over 5 • 10 4 
independent scenarios of users’ drops and channel coefficients. 
As seen, Vf approaches 1 for realistic values of P up to 
25% of the maximum rate, whereas for R = 30% the typical 
transmit power levels which are used in the uplink of present 
cellular systems, are not enough to ensure a high Vf. 

Next, we analyze the system performance in terms of 
EE and data rate. Fig. 2 shows the average GEE for the 
following resource allocation policies: (a) Algorithm 1 with 
R = 20%; (b) Algorithm 1 without QoS constraints, i.e. 
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Table II 

K = 5; M = 50; r = 10 -2 . AVERAGE NUMBER OF REQUIRED ITERATIONS TO REACH CONVERGENCE VERSUS Pmax FOR: (A) ALGORITHM 1 FOR GEE 
MAXIMIZATION WITH R = 20%; (B) ALGORITHM 1 FOR GEE MAXIMIZATION WITH R = 0%; (C) ALGORITHM 2 FOR DISTRIBUTED RESOURCE 
ALLOCATION WITH R = 20%; (D) ALGORITHM 2 FOR DISTRIBUTED RESOURCE ALLOCATION WITH R = 0%. 



Maximum power P [dBW] 


P = -38 

P = -34 

P = -30 

P = -26 

P = -22 

P = -18 

P = -14 

P = -10 

Algorithm 1. R = 0% 

2.63 

3.69 

4.68 

6.30 

6.53 

6.49 

6.50 

6.51 

Algorithm 1. R = 20% 

2.63 

3.67 

4.87 

6.68 

6.70 

6.76 

6.76 

6.77 

Algorithm 2. R = 0% 

1 

1.01 

1.07 

1.42 

2.54 

3.66 

4.12 

4.50 

Algorithm 2. R = 20% 

1 

1.01 

1.07 

1.42 

2.54 

3.67 

4.19 

6.71 



P [dBW] 


Figure 2. K = 5; M = 50; r = 1 0 “ 2 . Average achieved GEE versus P 
for: (a) Algorithm 1 for GEE maximization with R = 20%; (b) Algorithm 1 
for GEE maximization with R = 0%; (c) sum-rate maximization by adapted 
Algorithm 1; (d) Maximum transmit power allocation. 

R = 0%; (c) Algorithm 1 specialized to maximize the sum- 
rate; (d) Maximum power allocation, i.e. Pk = P for all k, 
considered as a baseline scheme. In scheme (a), if the problem 
turns out to be unfeasible, the QoS constraint is relaxed and 
the solution from scheme (b) is taken. For low values of P, 
this circumstance is very frequent, and indeed schemes (a) and 
(b) perform similarly. At the same time, schemes (b) and (c) 
also perform similarly for low P, thus suggesting that in this 
range of P, GEE and sum-rate maximization are equivalent. 
Instead, different performance are achieved for larger values 
of P. Indeed, increasing P eventually allows attaining the 
peak of the GEE. At this point, the GEE achieved by scheme 
(b) remains constant, as using the excess transmit power 
would only decrease the GEE. Indeed, the GEE achieved 
by scheme (c) decreases, because this scheme makes use of 
the excess transmit power to maximize the sum-rate. Instead, 
scheme (a) strikes a balance between these two extremes. 
Some of the excess power is used to fulfill the QoS constraints, 
which results in a slightly lower GEE. However, once the 
constraints are met, the transmit power is not further increased 
and the achieved GEE keeps constant. We remark that this 
slight reduction of the GEE grants a higher minimum rate. 
In particular, in the saturation region of the GEE, the average 
minimum rate granted by scheme (b) is 1.6 [bit/s/Hz/user], 



P [dBW] 


Figure 3. K = 5; M = 50; r = 10 2 . Average achieved minimum users’ 
energy efficiency versus P for: (a) Algorithm 1 for Min-EE maximization 
with R = 20%; (b) Algorithm 1 for Min-EE maximization with R = 0%; 
(c) Min-rate maximization by adapted Algorithm 1. 

whereas it increases to 2.35 [bit/s/Hz/user] when scheme (a) 
is used. 

Similar considerations can also be made when Algorithm 1 
is used to maximize the minimum of the users’ EEs. Fig. 
3 compares the minimum EE (with Wk = 1 for all k) 
versus P, achieved by: (a) Algorithm 1 with R = 20%; (b) 
Algorithm 1 with R = 0%; (c) min-rate maximization by 
adapted Algorithm 1. The results show a similar behavior as 
for Fig. 2. 

Fig. 4 compares the GEE performance of the centralized 
Algorithm 1 and of its distributed counterpart Algorithm 2 for 
R = 0% and 20%. We observe that, while the centralized 
scheme suffers a little performance gap in terms of GEE 
when QoS constraints are introduced, having minimum rate re¬ 
quirements causes a larger GEE degradation in the distributed 
scenario, especially for increasing P. This is expected because 
unlike the centralized scheme, in the distributed setting the 
interference among the users is not jointly managed, which 
results in high multi-user interference, especially for large P. 

Table II shows the average number of iterations required by 
Algorithms 1 and 2 to converge when R = 0% and R = 20%. 
Convergence is declared when ||qM — q( l - 1 )|| 2 /||q(*)|| 2 < 
Hr 1 . It is seen that both algorithms converge after a small 
number of iterations, which slightly increases for larger values 
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Figure 4. K = 5; M = 50; r = 10 ~ 2 . Average achieved GEE versus P for: 
(a) Algorithm 1 for GEE maximization with R = 20%; (b) Algorithm 1 for 
GEE maximization with R = 0%; (c) Algorithm 2 for distributed resource 
allocation with R = 20%; (d) Algorithm 2 for distributed resource allocation 
with R = 0%. 


Figure 5. K = 3; N = 16, M = 3. Average achieved GEE versus P for: 
(a) Algorithm 3 for GEE maximization with R = 20%; (b) Algorithm 3 for 
GEE maximization with R = 0%; (c) Algorithm 4 for distributed resource 
allocation with R = 20%; (d) Algorithm 4 for distributed resource allocation 
with R = 0%; (e) sum-rate maximization by adapted Algorithm 3 


of P. This is because increasing P results in a larger feasible 
set. Observe that the distributed algorithm exhibits faster con¬ 
vergence than the centralized one. This makes it particularly 
suitable for self-organizing networks. 


B. Relay-assisted OFDMA interference network 

Consider a relay-assisted multi-cell network as described in 
Section II-C2. Assume S = 3, N = 16, K = 3, M = 3, and 
B = 180 kHz. The UEs are placed at a distance from the relay 
randomly generated in the interval [100; 300] m and the same 
path-loss model and noise parameters as for previous figures 
are used. 

Fig. 5 reports the GEE versus P achieved by the following 
resource allocation algorithms: a) Algorithm 3 with R = 20%; 
b) Algorithm 3 with R = 0%; c) Algorithm 4 with R = 20%; 
d) Algorithm 4 with R = 0%; e) Algorithm 3 specialized to 
maximize the sum-rate. If the problems with QoS constraints 
happen to be infeasible, the solution of the corresponding 
scheme with R = 0% is taken. Similar remarks can be made 
as for the massive MIMO system. In particular, by introducing 
QoS constraints we trade-off a slight reduction in GEE with 
a significant increase of the users’ minimum rate, which, 
for the simulated scenario, increases from 3.8 [bit/s/Hz/user] 
when R = 0%, to 7.16 [bit/s/Hz/user] when R = 20%. As 
for the comparison between the centralized and decentralized 
approaches, we observe that the gap is rather limited in this 
scenario. This might be due to both the channel diversity 
granted by the use of multiple sub-carriers and to the rather 
limited number of non-cooperating transmitters. Indeed, it is 
known that that a limited number of players tend to result in 
a lower price-of-anarchy in game-theoretic settings [45]. 


C. Computational complexity discussion 

The computational complexity of the proposed algorithms 
depends on the number of iterations required to reach con¬ 
vergence as well as on the complexity of each iteration. 
A theoretical analysis for the number of iterations to reach 
convergence is very challenging. In fact, no theoretical results 
are available also for simpler scenarios than energy efficiency 
maximization. To overcome this problem, we have resorted to 
a numerical analysis. From Table II, it follows that the number 
of outer iterations for the single-resource-block case is very 
limited for both centralized and distributed algorithms (the 
latter require a slightly lower number of iterations). Similar 
results are obtained also in the multi-resource-block setting. 
The complexity of each iteration is analyzed in the following. 
The centralized algorithms are considered first. 

1) Centralized Algorithms: We begin by observing that 
the complexity of the centralized algorithms depends on the 
complexity of each fractional problem to solve. To fix ideas, 
let us consider the maximization of the GEE - similar consid¬ 
erations apply to the maximization of the Min-EE. Although 
there exist several different methods to solve fractional pro¬ 
gramming problems, the most widely used is the Dinkelbach’s 
algorithm. As shown in Algorithm 5, Dinkelbach’s algorithm 
works by finding the solution of an auxiliary problem in each 
iteration, say x*, and then updating the parameter A as A n = 
/(x*)/p(x*). One well-known result of the Dinkelbach’s 
algorithm is that the sequence {A,,},, converges with a super- 
linear rate. Moreover, this convergence rate does not depend on 
the complexity required to compute x*. In this respect, observe 
also that the auxiliary problems to obtain {x*} are convex and 
therefore can be solved with polynomial complexity in the 
number of variables and constraints (which are KN and 2K, 
respectively). Putting all these facts together, we can conclude 
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that the centralized algorithms have a polynomial complexity 
in each stage. This makes them of affordable complexity. 

2) Decentralized Algorithms: As for the distributed case, 
a similar analysis can be performed since the proposed al¬ 
gorithms consist of an outer loop, which requires to solve 
at each stage K fractional problems with N variables and 2 
constraints. Following the above arguments, we can state that 
the algorithms have a polynomial complexity in each stage. 
Observe that the complexity of each stage does not depend 
on K (as in the centralized case) but only on the number of 
resource blocks N and on the number of constraints (which 
is two). The parameter K only scales the complexity of each 
iteration in a multiplicative way. 

VIII. Conclusion 

In this paper, we have proposed a framework to develop 
centralized and decentralized power control algorithms for 
EE optimization in wireless networks. Unlike most previous 
related works, we have considered rate constraints and a more 
general SINR expression so as to encompass emerging 5G 
technologies. The resulting optimization problems have been 
tackled by an interplay of fractional programming, sequential 
convex optimization, and game theory. This has allowed us 
to derive centralized algorithms achieving first-order optimal 
points of the GEE and of the minimum of the users’ EEs, 
and to develop decentralized algorithms that have been shown 
to converge to the GNE of the associated game. The analysis 
above has been performed in the case in which a single or 
multiple resource blocks are used for transmission. Numerical 
results have been used to corroborate the theoretical results. 
For this purpose, two case-studies have been considered: a 
hardware-impaired massive MIMO network and a multi-cell 
multi-carrier relay-assisted network. 

The numerical analysis indicates that the centralized algo¬ 
rithms are quite robust to the enforcement of demanding rate 
constraints. Instead, the distributed algorithms are more sen¬ 
sitive to rate constraints, especially for increasing maximum 
feasible powers, due to the lack of centralized interference 
management. Also, the centralized algorithms perform better 
than their distributed counterparts, both with and without 
rate constraints, at the expense of a higher computational 
complexity and feedback requirements. 

An important follow-up of this work is to analyze the 
gap between the proposed centralized framework and the 
global optimum of the GEE and minimum EE. To this end, 
an optimization framework is required, which allows one to 
effectively determine the global solution of energy-efficient 
optimization problems in interference-limited networks. This 
is a challenging problem, which is still much open. An attempt 
in this direction is taken in [46], wherein tools borrowed from 
monotonic optimization theory, fractional programming, and 
sequential optimization are combined together to characterize 
the global maximum of the energy efficiency in interference- 
limited wireless networks. 

Appendix A 

Major concepts from fractional programming are reviewed 
here. For a more comprehensive overview, we refer to [16]. 


Definition 1 (Fractional program). Let C C R M be a convex 
set, and consider the functions f : C —> ffi. and g : C —> M + . 
A fractional program is the optimization problem 


max 

xec 


/(*) 


( 88 ) 


Proposition 8 ([39]). An x* £ C solves (88) if and only if 
x* = argmax xe c {/(x) — A*< 7 (x)}, with A* being the unique 
zero of F{ A) = max x6C {/(x) - Xg(x)}. 


This result allows us to solve (88) by finding the zero of 
F( A). An efficient algorithm to do so is the Dinkelbach’s 
algorithm [39], reported in Algorithm 5 for the reader’s conve¬ 
nience. If ,/(x) and g(x) are concave and convex, respectively, 
then the Dinkelbach’s algorithm requires to solve one convex 
problem in each iteration. 9 Moreover, the convergence rate of 
Dinkelbach’s algorithm is known to be super-linear [39], 


Algorithm 5 Dinkelbach’s algorithm 

Set e > 0; A = 0; 

repeat 

x* = arg max x gc {/(x) - Xg(x)} 
F = /( x *) - Xg(x*); 

A = /(x*)/fl(x*); 

until F < e 


A considerable extension of (88) is to consider the maxi¬ 
mization of the minimum of a set of ratios {/^(x)/< 7 i(x)}( =1 . 
This problem is usually referred to as generalized fractional 
programming, and has been first studied in [40], wherein an 
optimization procedure is provided, based on a modification 
of Dinkelbach’s algorithm. Specifically, the auxiliary function 
to be considered is F( A) = mini<j</ {/,(x) — A^i(x)}, and 
the algorithm works as shown in Algorithm 6. 


Algorithm 6 Generalized Dinkelbach’s 

algorithm 

Set e > 0; A = 0; 



repeat 



x* = argmaxx 

gc {mim<i<j [ffx ) 

1 - Afiff(x)]}; 

F = mini<j<j 

{/i(x*) - A5i(x*)}; 


A = mini<j</ 

/»(x*)/s*(x*); 


until F < e 




Then, if each ratio has a concave numerator and a convex 
denominator, we can solve the generalized fractional problem 
by solving a sequence of convex problems, 10 with a linear 
convergence rate [40]. 


Appendix B 

In the i-th iteration of Algorithm 1, we solve (41) and 
compute the corresponding vector pW of transmit powers, 
which maximizes the lower-bound at the i-th iteration, as 
given by (39), subject to the same constraints of the original 
problem (28). Then, the following chain of inequalities holds 

Ip( P (i) ) > <Mp W ) > = VKp^) (89) 

9 It is also required that A > 0 in each iteration. This can be shown to always 
hold if the algorithm starts with A = 0, and provided max® /(*) > 0. 

10 Observe that the minimum of concave functions is also concave [38]. 
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wherein the first inequality follows because ip, is a lower- 
bound of ip, the second inequality follows because pM is the 
maximizer of while the final equality holds because the 
parameters (i k > and b k 1 in ip, are such that the bound is tight 
in p(' i_ 1 i. As a consequence of (89), the value of ip increases 
after each iteration, and the algorithm must converge because 
ip is upper-bounded. 

Next, let us denote by ip and p the lower-bound of ip and the 
power vector at convergence. By construction, p maximizes 
ip subject to the same constraints of the original problem 

(28) , and therefore fulfills the associated KKT conditions. 
Such KKT conditions are the same as the KKT conditions 
of the original problem (28), except for the different objective 
function. However, upon convergence of Algorithm 1, we have 
ip( p) = ip{ p) and ViP(p) = Vippp), and the thesis follows. 

Appendix C 

The first part of the proof follows along the same lines of 
Proposition 1. As for the KKT conditions, we remark that in 
this case we can not directly consider the KKT conditions of 

(29) , because the objective of (29) is not differentiable. We 
can remove the non-differentiability by expressing (29) in its 
equivalent epigraph form: 

max i subject to ri k > t V/c. (90) 

{te R+.pe’P} 

Now, let us consider a modified version of Algorithm 1, which 
operates on (90) rather than on (29). In each iteration of this 
modified Algorithm 1, (37) is again used to lower-bound r) k 
with ijk■ This amounts to solving an approximation of (90), 
with fjk in place of i] k , and updating the parameters a!*' 1 
and b k 1 as in Algorithm 1. Following similar arguments as 
in the proof of Proposition 1, it follows that this modified 
version of Algorithm 1 converges to a point fulfilling the 
KKT conditions of (90), the only difference being that now the 
lower-bound is computed for the constraint function ij k rather 
than for the objective function. However, since rjk > fjk, the 
solution of the approximate problem is always in the feasible 
set of (90). Moreover, upon convergence, r) k and fjk are equal, 
and so are their gradients. Therefore, the modified version 
of Algorithm 1 yields a power vector fulfilling the KKT 
conditions of (90). Observe now that replacing ijk with fjk in 
(90) yields the epigraph form of (45), i.e., the problem which 
is solved in each iteration of Algorithm 1 by the Generalized 
Dinkelbach’s method. Hence, Algorithm 1 and its modified 
version (considered in this proof) converge to the same power 
allocation vector p, and hence the thesis. 


Appendix D 

Using [ 8 , Sect. II.A], it follows that solving the EE problem 
(47) is equivalent to finding the root of the nonlinear function 


$(A k ) = max {Slog (1 + 7 fe ) - X k {p c , k + p k )} (91) 

Pk 6R+ 


where X k £ R+ is such that $(A{.) = 0 and can be 
obtained through the Dinkelbach’s algorithm. Setting to zero 
the derivative of (91) with respect to p k yields: 


B d'y k 
1 + 7 k dpk 


— A k = 0. 


(92) 


From (15), using (14) one gets dik/dp k — Ffc(l~ =r ) 2 , from 
which (92) reduces to: 


Bp-k L _ 7fc\ 2 

1 + 7 k\ 7 k ) 


Xk = 0 . 


(93) 


In the attempt of solving (91), let us study the properties 
of (93) as a function of 7 k . This amounts to studying 
fkilk) = 0 with fkilk) = - (2B/i fe 7 fc + X k ^ 2 k ) i k + 

7 j" (B/i k — X k ). The first derivative of fkilk) with respect to 

7 k is 

dJ ^ = 2 Bp k ilk - Ik) - Afc7fc < 0 (94) 


where the last inequality follows since 0 < 7 /,. < j k , and 
Hk > 0, A k > 0, and 77 . > 0. As a consequence, fkilk) is 
a strictly decreasing function, which admits a solution if and 
only if fki 0 ) > 0 and fkilk) < 0 (with the equalities not 
simultaneously active). Since fkilk) = ~^kl k (1 + 7fc) < 0, 
we need to ensure fki 0) = i\ (B/r^ — Xk) > 0. This trans¬ 
lates into Bfik > Xk- By solving the second-order equation 
fkilk) = 0, one gets i k = v k (A fc ) where u k (x) is defined in 
(58). Plugging 7 k = v k (Afc) into (14) yields (60). 


Appendix E 

Rewrite p k ,n in (11) as n k , n = Oik, n /i^l,n + h,n) where 
Ik,n = J2jjtk u kj,nPj,n- Using [17, Theorem 4], the GNE is 
unique if ||<9Ifc/<9p_fc|| sup Ifc6R iv ||Bfc(p_fc)/dIfc|| < 1, for all 
k. with Ifc = [7fc ) i,... ,/fc,Ar] T - From [17, Eq. (19)], we have 


||<9Ifc/dp_fc|| 


I< N 


\ n =1 


(95) 


|| <9£fc(p_fc)/<9Ifc|| 


/ N N 

I E E 

1=1 n—l 


\dn k ,niX' k ) /dlk,e\ 2 - 


(96) 


Observe now that n k , n iX' k ) > 0 when ji kn > X' k (see 
Appendix D). After lengthy computations (not shown for space 
limitations), one gets 


9p*k, n 

dlk,e 


^fc,nlin=r 


1k,r 


9k,n (A k ) 


fk,t 




(97) 


where we have defined (for notational compactness) s k , n as 
in (83), whereas £ k ,e is defined in the text of Proposition 7. 
Plugging (97) into the RHS of (96) yields 


N 


ees* nes* 


T k,n 

9k,n (A fc ) 


(98) 


where S k = {n = 1 ,..N : fi k , n > X' k }. Putting all the 
above results together. Proposition 7 follows. 
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